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We have computed the low energy quantum states and low frequency dynamical susceptibility of 
complex quantum spin systems in the limit of strong interactions, obtaining exact results for system 
sizes enormously larger than accessible previously. The ground state is a complex superposition of 
a substantial fraction of all the classical ground states, and yet the dynamical susceptibility exhibits 
sharp resonances reminiscent of the behavior of single spins. These results show that strongly 
interacting quantum systems can organize to generate coherent excitations and shed light on recent 
experiments demonstrating that coherent excitations are present in a disordered spin liquid £jj . The 
dependence of the energy spectra on system size differs qualitatively from that of the energy spectra 
of random undirected bipartite graphs with similar statistics, implying that strong interactions are 
giving rise to these unusual spectral properties. 



The motivation for understanding the dynamics of 
quantum systems continues to grow as technology minia- 
turizes. That quantum and classical dynamics differ fun- 
damentally is underscored by the indications that quan- 
tum dynamics can be harnessed in the form of quantum 
computation to solve problems that are intractable using 
the processes of classical physics Q. 

The promise and difficulty of quantum dynamics both 
arise because the amount of information needed to spec- 
ify a quantum state grows exponentially with the system 
size. The number of complex numbers needed to spec- 
ify a general quantum state of N quantum-mechanical 
two-state spins (or qubits) is 2^. A classical computer 
performing 10 floating point operations per second |3( 
would take 10 8 years to perform a single operation when 
N = 100 and 10 21 years when N = 144. 

Here we study a frustrated spin model that is a quan- 
tum generalization of the two-dimensional ± J Edwards- 
Anderson (E- A) spin glass model Q , a canonical example 
of a classical system whose competing interactions give 
rise to many low-energy states. The essential physical in- 
gredients of the E-A model arise in a wide variety of opti- 
mization problems in many fields 5] : the system cannot 
satisfy simultaneously all its constraints, and many dif- 
ferent configurations are equally effective in minimizing 
the energy. Though the two-dimensional ± J E-A model 
is simpler than the three-dimensional version — it is dis- 
ordered at all nonzero temperatures [y] and individual 
ground states can be found in a time that scales as a poly- 
nomial of the system size Q — it has a large ground state 
degeneracy and a complex energy landscape H, H, @- 
Motivated by recent experiments on LiY x Hoi_ x F4 that 
demonstrate that quantum tunneling has profound ef- 
fects on the dynamics of a three-dimensional Ising spin 
glass with dipolar couplings 0, 0, Q], we study the 



quantum system obtained by adding a small quantum 
tunneling term to the two-dimensional E-A model. 

We calculate numerically exact eigenstates and dynam- 
ical response functions of quantum E-A models with up 
to 144 spins in the limit of strong interaction strength. 
At low frequencies the excitation spectrum is remark- 
ably sparse. The number of excitations per unit energy 
at low energies does not vary significantly with system 
size, even though the number of eigenvalues grows expo- 
nentially with system size while the energy bandwidth is 
only growing polynomially. The eigenvalues and eigen- 
vectors of each realization are obtained by diagonalizing 
the adjacency matrix of a graph, and we demonstrate 
that strong correlations play a vital role by comparing 
the spin glass excitation spectra to spectra of graphs 
with similar connectivity statistics but randomly chosen 
connections. These theoretical results provide a natural 
framework for understanding recent experiments demon- 
strating the presence of sharp, saturable resonances in 
the quantum spin liquid LiYo.955Hoo.o45F4 0. 

We study two-dimensional systems with periodic 
boundary conditions in which spin-1/2 spins interact 
with nearest neighbors on V7V x V^V square lattices. The 
quantum Hamiltonian is 
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The sum (ij) is over all nearest 



neighbor pairs. A given sample has a fixed realization 
of bonds in which each bond Jij is chosen to be — J and 
+ J with equal probability. This Hamiltonian, which has 
been studied by many groups |l2T | , is believed to be rele- 
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vant to the experimental system LiY x Hoi_ x F4 [H HclllH . 
where the coefficient T is tunable because it is propor- 
tional to an applied transverse magnetic field. The r — > 
limit is the two-dimensional ± J Edwards- Anderson (E- 
A) model 0]. 

As r/J — > the low energy quantum states \i/j n ) be- 
come superpositions of states corresponding to ground 
state configurations \a) of the classical model, 



(2) 



The number of ground states of the classical two- 
dimensional ±J E-A model grows with N much more 
slowly than the number of configurations (though still 
exponentially) 0, 0] ; we use an algorithm that we have 
developed that finds all the classical ground states effi- 
ciently and standard degenerate perturbation the- 
ory [l4| to compute the low energy quantum states of 
about 90% of 10 x 10 realizations and about 40% of 
12 x 12 realizations 0. In the limit r/J — > 0, each 
low energy quantum eigenstate is a superposition of clas- 
sical ground state configurations related to each other by 
serial flipping of individual flippable spins, where a flip- 
pable spin is one with an equal number of satisfied and 
unsatisfied bonds ^(|. This follows because the matrix 
element (a\ a ix) 1/3) of the quantum tunneling term 
between any two classical configurations corresponding to 
the quantum basis states \a) and \(3) is nonzero only if the 
two states differ by a single spin flip, and, to lowest order, 
only states \a) and \(3) with the same energy contribute. 
The number of flippable spins clearly is no greater than 
N, the number of spins in the system, so the matrix 
characterizing the possible transitions between classical 
ground states is extremely sparse, and thus well-suited 
for diagonalization using Lanczos techniques \vH Hil . 
We have computed the dynamical magnetic susceptibility 
by finding low-energy eigenvalues and eigenstates using 
the sparse matrix ARPACK numerical library [^j with 
C++ bindings [2]]]. Most of the results shown here are 
for individual realizations; though there are large samplc- 
to-sample variations in the number of ground states for 
a given system size the qualitative results on which 
we focus are robust. 

The ground state dynamical magnetic susceptibility 
X {oj) characterizes the response of a system at zero tem- 
perature to a magnetic field applied along the z axis os- 
cillating at frequency oj .22]. The susceptibility consists 
of sets of Dirac 6- function peaks (that in physical systems 
spread out into a finite width in frequency due to deco- 
herence processes) ; each peak occurs at a frequency that 
is H times the energy difference between an excited state 
and the ground state. For the ± J spin glass, as T / J — > 
the value of J affects only the energy zero and the sus- 
ceptibility at frequencies oj satisfying hui <C J depends 
only on the ratio ftw/Y. 

Figure^shows the zero temperature dynamic magnetic 
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FIG. 1: Zero temperature dynamic magnetic susceptibility 
X (oj) of systems of size 6 x 6, 8 x 8, 10 x 10 and 12 x 12. 
The peaks in the susceptibility occur at frequencies 0) that 
satisfy hcu — E n — Eo, where E n is the energy of an excited 
state and Eo is the energy of the ground state. The density 
of low-energy excitations does not increase appreciably as the 
system size increases, even though an exponentially increas- 
ing number of states are in an energy bandwidth that grows 
approximately linearly with the system size. 



susceptibility of systems of size 6x6,8x8, 10 x 10, and 
12 x 12. The density of low energy excitations increases 
extremely slowly with system size. This result is sur- 
prising because the number of energy eigenvalues grows 
exponentially with the number of spins N, while these 
energies all lie within a bandwidth that grows roughly 
linearly with N. 

To obtain context for these results, we interpret the 
Hamiltonian matrix for the quantum spin glass as the 
adjacency matrix of an undirected bipartite graph [2^ . 
24, 25] in which each classical ground state is a node and 
edges connect every pair of classical ground states cou- 
pled by the quantum term in the Hamiltonian. The graph 
is bipartite because the edges connect states that differ 
by a single spin reversal, one of which has an even and the 
other an odd number of up spins. The spin glass graphs 
have a modest number of disconnected pieces, called clus- 
ters [2(| . Figure |21 compares the density of energy levels 
of the largest cluster of a 10 x 10 spin glass realization 
(with 17040 nodes and 77684 edges) to the density of en- 
ergy levels of a symmetric bipartite random matrix with 
10000 nodes and 50000 edges. The bipartite random ma- 
trix has a large energy gap between the ground state and 
first excited state, and once this gap is exceeded the den- 
sity of energy levels is much greater than at low energies 
in the spin glass. The energy level spacing between the 
excited states of the random matrix is approximately in- 
versely proportional to the number of nodes. 

We have compared other properties of the graphs 
underlying the quantum spin glass to those of random 
bipartite graphs. The degree distribution [2^ describ- 
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FIG. 2: Density of states as a function of energy at low ener- 
gies for the largest connected component (with 17040 nodes 
and 77684 links) of the graph characterizing a 10 x 10 spin 
glass realization and of a bipartite random matrix with 10000 
nodes and 50000 links. The quantum ground state energy Eq 
for the spin glass is Eo/Y = —10.1949, and for the bipartite 
graph Eo/Y = —10.2335. The ordinate shows the number of 
eigenvalues in a bin of width 0.01. The bipartite random ma- 
trix has a large gap between the ground state and first excited 
state, and, once the gap is exceeded, a much larger density of 
states than the spin glass. 



ing the number of links emanating from the nodes of 
the spin glass graphs is even narrower than the Pois- 
son distribution of a bipartite random graph with the 
same mean degree. Figure [3] shows the clustering coeffi- 
cient C [2^, which for bipartite graphs is the probabil- 
ity that two nodes with a common second neighbor are 
themselves second neighbors 30] . The clustering coeffi- 
cients of spin glass graphs are significantly larger than 
those of bipartite random graphs with the same num- 
ber of nodes and edges |30j , and are close to those of 
graphs describing N noninteracting spins, which have 2 N 
nodes, each node with degree N, and clustering coeffi- 
cients C = 4/(AT+ 

Though some statistical properties of the spin glass 
graphs are similar to those of graphs for noninteract- 
ing quantum spins, the ground state of the quantum 
spin glass differs significantly from that of noninteracting 
quantum spins. Figure^] (inset) displays the spins in the 
quantum ground state of a 12 x 12 system that are fluctu- 
ating, in that |(<Si Z )| 7^ 1, where (5j Z ) is the expectation 
value of the zth component of the i th spin. The fluctu- 
ating spins form connected "bunches" |l3j with up to 15 
spins, so a description in terms of noninteracting spins 
is not appropriate. The main panel of Figure 01 shows 
the participation ratios P a of different classical ground 
states in the quantum ground state, where P a — \c a o\ 2 , 
with the c a Q defined in Eq. The spin glass partici- 
pation ratios vary over several orders of magnitude (note 
the logarithmic abscissa), while for a system of N non- 
interacting spins, the participation ratio is 1/2 N for all 
classical configurations. Figure 0] shows that the partici- 



FIG. 3: Clustering coefficients C of spin glass graphs, of 
bipartite random graphs with the same number of nodes and 
edges, and of graphs for noninteracting quantum spins, versus 
number of nodes. The clustering coefficients of the spin glass 
graphs are significantly larger than those of random bipartite 
graphs, and close to those of graphs for noninteracting spins. 
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FIG. 4: Histogram of the probability distribution of classical 
ground states whose participation ratio P a is in a given range 
in the quantum ground state, on a logarithmic scale. Here, 
Pa = |ccto| 2 , where the c an are defined in Eq. J2J. The distri- 
butions of participation ratios differ significantly between the 
spin glass and the random bipartite matrix, though in both 
cases the P a vary over many orders of magnitude (note the 
logarithmic abscissa). Inset: Spins denoted with open cir- 
cles have mean magnetization {Si) satisfying |(iS»)| 7^ 1 in the 
quantum ground state. 



pation ratios of the adjacency matrix of a bipartite ran- 
dom graph also vary over several orders of magnitude, 
but with quite different statistics. 

We expect the perturbative methods we use to be valid 
so long as the energy arising from the quantum perturba- 
tion TN) is smaller than the energy gap between the 
classical ground state and lowest classical excited states 
(= J), so that the procedure is valid only for V ^ J/N. 
Exact diagonalizations of very small (3x3 and 4x4) sys- 
tems are consistent with this expectation. 

Our results, particularly the low density of well-defined 
low-frequency excitations, provide a framework for un- 



4 



derstanding the generation of coherent, saturable excita- 
tions in a complex quantum spin liquid 0. Saturability 
arises naturally when the density of excitations is low 
because a frequency that induces a transition from the 
ground state to an excited state will not be able to induce 
a second transition from the excited state. When exter- 
nal excitations couple only two states, the system will 
display coherent oscillations and will have the capability 
of encoding phase-coherent information 0]. An impor- 
tant future goal is to understand how to implement con- 
trolled dynamics involving more than two energy levels 
in this strongly interacting system, as has been achieved 
in the weakly interacting regime in the context of NMR 
quantum computation |32j |. 
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